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Improved theoretical predictions for tlie fine-structure splitting of 2^P./ levels in helium are obtained by the 
calculation of contributions of order ot' Ry. New results for transition frequencies, i^oi = 29 616 943.01(17) 
kHz, and i/yz — 2 291 161.13(30) kHz, disagree significantly with the experimental values, indicating an out- 
standing problem in bound state QED. 
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The fine-structure splitting of the helium states is an 
intrinsically relativistic effect, and arises from the interaction 
of spins and orbital angular momentum. The value of this 
splitting has been measured with increasing precision over the 
last years 01 1^ 1^1 IS 1^ . Since the splitting is proportional 
to Ry, these accurate measurements make helium a can- 
didate for determining the fine structure constant a, provided 
that the higher order in a corrections can be sufficiently well 
understood. The most accurate determination of a at present 
comes from the 5 — 2 of the electron. This determination de- 
pends sensitively on complicated multi-loop calculations per- 
formed by Kinoshita and by Remiddi and coworkers |7], and 
therefore requires independent confirmation. In response to 
significant experimental effort 01 S S 01' we present 
here the calculation of the Ry contribution to helium fine 
structure, so that these experiments can be used to provide an 
independent determination of a. 

Several recent advances in bound state Quantum Electro- 
dynamics (QED) have made the calculation of higher order 
corrections to helium fine structure possible. Specifically, 
Yelkhovsky in Ref. 1 8] has shown how to use dimensional reg- 
ularization in the calculation of helium energy levels, and to- 
gether with Korobov has obtained in |9] numerical values for 
the Ry contributions to the ground state. Next, in Ref. tlOll 
a Foldy-Wouthuysen transformed QED Langrangian was used 
to derive all effective Ry operators for arbitrary states of 
few electron atoms. More recently, together with Jentschura 
and Czarnecki, we have obtained in Ref. Ill 111 general for- 
mulae for Ry coiTection to hydrogenic energy levels, in- 
cluding the fine structure. The calculational approach of these 
works 1 9, 11] and the present paper is based on dimension- 
ally regularized QED. The parameter e, related to the space 
dimension = 3 — 2 e, plays the role of both infrared and 
ultraviolet regulator, as some Ry terms are divergent in 
d = 3 space. This artificial parameter e is used to derive var- 
ious terms, and we will explicitly demonstrate its cancella- 
tion in their sum. Natural relativistic units will be used with 
?i = c = eo = m = 1, so that = 4 tt a. 

The fine structure in order m (a^ Ry) can be written as 
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where is the Bethe logarithmic correction of Eq. ( I15> , 



and is an effective Hamiltonian of order ma*. We will 
concentrate in this work on a complete derivation of i/*^^), as 
the other terms contributing to order m {ot' Ry), E^ and 
the second order term called Eg, have already been obtained 
in 1 13]. Important terms of order moi' In a first calculated 
in Ref. I15II are confirmed in the present calculation. H'^'^^ 
consists of exchange terms and radiative corrections, where 
a photon is emitted and absorbed by the same particle. We 
consider first the exchange terms. Their derivation in general 
is quite complicated. We note that only two-photon exchange 
diagrams contribute and there are no three-body terms, which 
is a result of an internal cancellation. A feature of the calcu- 
lation that leads to considerable simplification is the fact that 
the order being calculated in nonanaly tic in . For example, 
i/*^^-' includes only two terms 

i?''' = — + — ^ [6Hri) + '5'(r2)] , (2) 
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and they can be derived from the two-photon exchange scat- 
tering amplitude. Similar results hold for the spin dependent 
m terms. If represents an effective Hamiltonian, it 
has to give the same scattering amplitude as in full QED. 
Therefore, we obtain the exchange contribution SH from the 
spin dependent part of the two-photon scattering amplitude, 
which is 
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where q = p'l — pi = P2 — P2- If one expands this amplitude 
in small external momenta one obtains 
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where u'-' 
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divergences cancel out with the low energy part where photon 
momenta are of the order of the binding energy. This low 
energy contribution gives the Bethe logarithm, described later 
in Eq. (I15> . and the correction 
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which is the transition term from dimensional regularization 
to the direct A — m{Z a)^ A cut-off in the photon momenta. 
Here 5 denotes the first order correction to iJ and E due to 
the spin dependent part of the Breit-Pauli Hamiltonian ij'*'. 
The resulting correction is a sum of two terms. The first one 
contributes to (7J('*V(-Eo - Hq)' H'^^^) in Eq. Q, and the 
second term is the effective Hamiltonian 



52H 



\n[{Zay 



[i ai{p[,pi) 



"5 1 2 

9 37 ^ 3 
+icr2{p'2,P2) + 2i(Ji{p2,P2) + 2i(T2(pi,pi) 
-cri{3,q)<^2{j,q)], (6) 

where we omitted a In 2A term. Together with Eq. (0} this 
gives the complete contribution due to exchange terms. When 
calculating expectation values on ^Pj states further simplifi- 
cations can be performed. Namely, the expectation value of a 
Dirac delta function with both momenta on the right or on the 
left hand side vanishes. Moreover, the nonrelativistic wave 
function is a product of a symmetric spin and an antisymmet- 
ric spatial function. This means that the expectation value of 
(Ti is equal to that of (T2- As a result the total exchange contri- 
bution He ~ 5iH + S2H is 
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The treatment of the radiative correction is different. We 
argue that radiative corrections can be incorporated by the use 
of electromagnetic formfactors and a Uehling correction to the 
Coulomb potential 
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The possible additional corrections are quadratic in electro- 
magnetic fields: see Ref. |11]. However, terms formed out 
of E, B, p, a can contribute only at higher order and thus can 
be neglected. Corrections due to the slope of formfactors and 
the vacuum polarization are obtained analogously to the Breit- 
Pauli Hamiltonian i/^"*), by modifying electromagnetic ver- 
tices and the photon propagator The result is 

= nZaiF[+2Fi + F{,)z[a,iplPi) + ^2{plP2)] 
-TT a{2 F[ + 2F!, + F^)i [ai {p[ ,pi) + a2 {p'2 , P2)] 
-2Tra{2Fi+F^ + Fi.)i[aiip'^,P2) + cj2 {p[ , Pi )] 
+7r a{2 F[ + 2 + di (j, q) ai {j, q) , (9) 

where by p" we denote momentum scattered off the Coulomb 
potential of a nucleus, and F' = F'{Q). There is also a low- 
energy contribution which is calculated in a way similar to this 
in Eq. (|5}, namely 
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The resulting effective Hamiltonian is 
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The complete radiative correction is a sum of Eqs. (|9} and 
(II It . namely Hn — S3H + S^H. Using symmetry 1 ^ 2 it 
takes the form 
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It is convenient to consider a sum of Eqs. Q and M2\ . as 
several logarithmic terms cancel out 
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The logarithmic terms agree with Refs. while non- 

logarithmic terms Qi are presented in Table I. 
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TABLE I: Operators due to exchange diagrams, slope of formfac- 
tors and the vacuum polarization, in atomic units with a prefactor 
ract j-K. The singular ^ dr/r integral is defined with an implicit 
lower cut-off A and the term In A + 7 is subtracted out. 

Operator t'oi[kHz] 1^12 [kHz] 



Qi = 


Zipi X <5^(ri)pi -ai 


2.854 


5.709 


Q2 = 




10.886 


-4.355 


Q3 = 


- ^ f ■ (Tl f • (T2 


4.132 


-1.653 


Qi = 


Sfjf ipi X(5^(r)pi-ai 


5.186 


10.372 


Qs = 




-1.328 


-2.656 






21.731 


7.418 



The remaining contribution is the anomalous magnetic mo- 
ment correction to the spin dependent operators. We derive it 
with the help of a nonrelativistic QED Hamiltonian obtained 
by a Foldy-Wouthuysen transformation of the Dirac Hamilto- 
nian including the magnetic moment anomaly n IIiHi 

-'2 -^4 
Hfw = y + e^*'-|(l + K)a-^-y 

-- {1 + 2k)['V ■ E + a ■ {E X n -TT X E)] 
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All the m operators obtained by Douglas and Kroll (DK) 
in JItIi can also be obtained from this Hamiltonian in Eq. ( I14t . 
see Ref.fl^. The anomalous magnetic moment operators are 
derived in a very similar way. They differ (see Table II) only 
by multiplicative factors from the DK operators. There is a 
one to one correspondence with Table I of Ref. 1 18] with 3 
exceptions. The operator Hs from Table II canceled out in 
DK calculation. The other two exceptions are related to the 
different spin structure of the next to last term in Eq. ( I14t . 
which leads to operators Hiq and Hu in our Table II. 

Apart from the Hi and Qi operators, second order con- 
tributions and low energy Bethe-logarithmic type corrections 
contribute to the fine structure. These contributions have al- 
ready been considered in our former work 1 13]. The second 
order contribution Es, beyond the anomalous magnetic mo- 
ment terms is the second term of Eq. Q, and the low energy 
contribution El is 
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where 5{. . .) denotes the correction to the matrix element 
(. . .) due to ij'"*^. Numerical results for all these contribu- 
tions is presented in Table III. 



TABLE II: Operators due to magnetic moment anomaly in atomic 
units with the prefactor m ol' /-k 

Operator j^oi [kHz] P12 [kHz] 



H2 -- 
Hz -- 
Hi -- 
H5 -- 
He -- 
Ht -- 
Hs -- 
Hj -- 
Hio 
H^i 

H\2 

Hn 
Hu 

His 
Hie 

H\7 



-¥|'x 7^-^i i^^-P^) 

3 Z f =! ri -f 



2 r 

3 

4r 



r X P2 • O"! 

r ■ 3\r ■ (72 



\V\^ X pi ■ CTl 

73- X P2 ■ CTl 

-;§7| xp2-ai 



2 

4 r- 



\p\^r-p2rxpx-Si 



r X (r ■ p2)pi ■ ui 

f X (r X pi ■ CTi) P2 • (?2 
pi • (?2 P2 • ffl 



8r 
1 

' 8r 



= j^Pi ^r- aif- a2 



ipl^{P-ff2 



Eh 



P2 ■ <?i + (f ■ CTi ) 

X (p2 ■ CT2) - ^ r ■ ^ir ■ a2f ■ P2) 
^ -jPi ■ SiPi X ■ P2 
= I pi • CTi (-pi ■ (T2 

+ 3pi • f p;- • CT2) 
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3.239 
0.267 
0.332 
0.749 
2.638 
-0.807 
-1.237 
-0.460 
0.093 
-0.376 
-0.193 
-0.447 
-14.908 
4.411 

4.618 
-0.483 

-1.643 



6.478 
0.534 
-0.133 
1.498 
-1.055 
-1.614 
-2.474 
-0.920 
0.187 
-0.752 
0.077 
0.179 
5.963 
-1.764 

-1.847 
-0.967 

0.657 



-4.208 



4.047 



TABLE III: Summary of contributions to helium fine structure, 
and iS'®-* including nuclear recoil corrections and the elec- 
tron anomalous magnetic moment at the level of the Breit- 
PauU Hamiltonian; = 137.03599911(46), m^/m^ = 

1.37093355575(61)10"*, Ry c = 3.289841960360(22) 10^^ Hz. 
Not indicated is the uncertainty due to a, which is 0.20 kHz for 1^01 • 
The last row includes the most recent experimental values. 





1^01 [kHz] 


ui2[kUz] 


Ref. 


Eq 


21.73 


7.42 




Eh 


-4.21 


4.05 




Es 


11.37(02) 


-1.25(01) 


im 


El 


-29.76(16) 


-12.51(27) 


[13] 




-0.87(16) 


-2.30(27) 






82.59 


-10.09 






-1557.50(06) 


-6544.32(12) 


[ISfiil, (49,201 


ew 


29 618 418.79(01) 


2 297 717.84 


[14],[19] 


total 


29 616 943.01(17) 


2 291 161.13(30) 




Drake 


29 616 946.42(18) 


2 291 154.62(31) 


|44] 


exp. 


29 616 951.66(70) 


2 291 175.59(51) 





Since all relevant contributions to helium fine structure 
splitting now seem to be known, we are at a position to present 
final theoretical predictions, which is done in Table III. Al- 
though we have included all terms up to order m a^, theoret- 
ical predictions are in apparent disagreement with the mea- 
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surements, as can be seen from the last row of Table lllll Let 
us analyze possible sources of this discrepancy. The numeri- 
cal calculation involves variational nonrelativistic wave func- 
tion. The parameter which controls its accuracy is the non- 
relativistic energy. Our wave function, consisting at maxi- 
mum of 1500 explicitly correlated exponential functions re- 
produces energy with 18 significant digits in agreement with 
the result of Drake in [14J. Matrix elements with this wave 
function are not as accurate as nonrelativistic energy, but they 
are sufficiently accurate for leading fine structure operators, 
and results agree with more accurate and independent calcu- 
lation of Drake in (3- For example, E^'*'^ agrees to 0.01 kHz 
and to 0.1 kHz. In fact almost all numerical calcula- 
tions have been performed by us and by Drake independently 
with one exception; we have not obtained recoil correction to 
the second order matrix element with Breit operators in E'-^'> . 
More important is the complexity of the derivation of m 
operators, namely Hi and Qi. We purposely derived Hi in a 
way very similar way to the derivation of the D-K operators to 
avoid accidental mistakes. We note that the Qi operators were 
obtained from the one-loop scattering amplitude in an almost 
automatic way, in contrast to the former very lengthy deriva- 
tion of Zhang I 2L.22..23J . with which we are in disagreement 
(see summary of Zhan g result s in Ref. (3)- In our previous 
papers with Sapirstein 1 13, 20] we pointed out several compu- 
tational mistakes and inconsistencies in Zhang's calculations, 
and therefore we consider the result of Drake, see Table [TO] 
to be incomplete. While it is possible that we have made 
a mistake somewhere, the other probable explanation of the 
discrepancy with experiments is the neglect of higher order 
terms, namely m a^. An indication of their importance is the 
recoil correction to the second order contribution, obtained 
by Drake in lll4ll . In spite of the small electron-alpha parti- 
cle mass ratio, this correction is very significant; for example 
61^01 = —10.81kHz. The mass ratio TOe/mc « 0.00014 is 
not much different from « 0.000053, and for this reason 
one can expect that iteration of Breit-Pauli Hamiltonian in the 
third order might also be significant. However, most m 
operators should be negligible, as i?*^^^ is akeady at the few 
kHz level, so that an additional power of a will make these 
operators contribute well below the experimental accuracy. 

In summary, we have obtained the complete Ry contri- 
bution to helium fine structure splitting. Theoretical predic- 
tions, including this result, are in disagreement with measure- 



ments fl','?,'?,'?,'?,'^. Therefore the determination of a from 
helium spectroscopy requires both checking the calculation of 
E^'^'> and the reliable estimation of the higher order i?^^) con- 
tribution, which is a challenging task. Therefore, at present, 
helium fine structure splitting is not competitive with respect 
to other determinations of a, for example, from the recent ex- 
periment on the photon recoil i24ll . 

I am grateful to Andrea Ferroglia and Michal Czakon for 
useful advice on evaluation of d-dimensional Feynman dia- 
grams. 
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